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I.   Introduction 

Present  day  platforms,  in  water  depths  from  10  to  1000 
feet,  are  primarily  of  welded  steel  jacket  construction.   As 
the  present  limit  of  about  1000  feet  was  approached  many 
problems  in  fabrication  and  installation  were  encountered, 
indicating  the  need  for  a  shift  to  a  new  technology.   More 
significantly,  with  inceasing  water  depth  the  natural  per- 
iods of  jacket  platforms  increase  from  about  one  second  in 
100  feet  of  water  to  approximately  four  seconds  in  1000  feet 
of  water.   There  is  significant  wave  energy  with  a  period  of 
four  seconds  and  dynamic  amplification  becomes  a  problem. 
The  expected  life  of  the  structure  is  shortened  due  to  high 
cycle  -  low  stress  fatigue  damage. 

One  alternative  structure  is  a  tension  leg  platform 
(TLP) .   A  TLP  (see  Figure  1-1)  is  a  bouyant  structure  held 
in  position  by  tension  members  secured  to  the  ocean  bottom. 
Because  the  platform  is  not  rigidly  connected  to  the  bottom 
its  natural  modes  of  vibration  are  significantly  different 
from  a  jacket  structure  of  similar  height.   A  TLP  in  a  water 
depth  on  the  order  of  1000  feet  will  have  natural  periods 
in  heave,  pitch,  and  roll  of  about  two  seconds,  well  out  of 
the  region  of  high  wave  energies. 


Since  the  TLP  has  only  recently  been  accepted  as  a 
potential  alternative  for  deep  water  drilling  there  are  few 
studies  of  its  dynamic  response  to  wave  forces.   Such  preli- 
minary studies  are  necessary  to  provide  information  for  the 
design  of  the  operating  platforms.   Response  prediction  re- 
quires specification  of  the  modal  exciting  forces,  for  example 
in  pitch,  roll,  and  heave.   For  low  everyday  sea  states  the 
exciting  forces  and  structure  responses  can  be  considered  to 
behave  in  a  linear  fashion.   Linearity  of  excitation  and  re- 
sponse suggests  that  a  frequency  domain  solution  technique 
is  appropriate;  especially  when  one  considers  that  ocean 
waves  are  a  random  process  most  often  described  by  a  wave 
amplitude  spectrum.   This  thesis  presents  an  estimate  of  the 
modal  exciting  force  spectra  in  heave,  pitch,  and  roll  in  a 
way  which  accounts  for  structure  geometry  and  wave  spreading. 


Figure  1-1  Tension  Leg  Platform 


II.   Problem  Statement 

The  modal  exciting  forces  on  a  tension  leg  platform  (TLP) 
will  be  dependent  upon  the  frequency  and  the  angles  of  inci- 
dence of  the  incoming  waves.   For  some  angles  of  incidence 
(at  a  given  frequency)  the  forces  on  the  various  legs  will 
act  together  to  produce  a  larger  total  force  or  moment.   For 
other  angles  of  incidence  the  phase  differences  will  cause 
the  forces  or  moments  exerted  on  the  various  legs  to  cancel, 
tending  to  produce  a  smaller  total  force. 

In  this  report  the  total  modal  force  of  the  platform  in 
heave,  pitch,  and  roll  will  be  divided  into  two  parts.   The 
first  will  be  the  force  exerted  upon  a  single,  isolated  leg 
by  waves  of  the  given  frequency.   The  second  will  be  a  geo- 
metric term  to  take  into  account  the  effect  of  the  phase 
differences  between  forces  on  each  of  the  legs.   This  second 
term  will  be  specifically  calculated  for  a  square  TLP  noting 
that  the  method  could  be  used  for  any  shape. 

The  calculations  so  far  imply  a  unidirectional  wave 
amplitude  spectrum.   Rarely  is  this  observed  in  practice. 
A  reasonable  wave  spectrum  can  be  approximated  as  being 
centered  about  one  bearing  from  the  structure  and 
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o 
diminishing  to  zero  at  bearings  90   either  side  of  the 

central  bearing.   Two  types  of  wave  spreading  functions 
will  be  discussed.   The  first  is  the  cosine  squared  spread- 
ing law  which,  although  widely  used,  allows  little  if  any 
leeway  in  varying  the  spreading  to  fit  actual  measurements. 
The  second  is  the  Gaussian  distribution  which  does  allow 
continuous  variation  of  spreading  simply  by  changing  the 
standard  deviation. 

The  modal  wave  force  spectra  upon  the  structure  will 
be  the  integral  over  all  incidence  angles  of  the  directional 
wave  amplitude  spectra,  weighted  appropriately  by  the  plat- 
form heave  force,  roll  or  pitch  moment  transfer  functions. 
The  resulting  spectra  will  be  normalized  in  this  thesis  by 
the  results  obtained  when  a  uniformly  spread  wave  spectrum 
is  used.   This  is  done  to  demonstrate  the  effect  of  wave 
spreading  on  the  heave,  pitch,  and  roll  exciting  force 
spectra. 
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III.  Formulation  of  Equations 

A.   Platform  transfer  function 

1.  Single  leg  transfer  function.   The  single  leg  wave 
amplitude  to  vertical  force  transfer  function,  re(co)/  is  a 
function  only  of  wave  frequency.   Because  the  leg  has  circu- 
lar cross  sections  there  is  no  dependence  on  0 ,  the  angle  of 
incidence.   This  thesis  assumes  that  rs(u>)  is  available, 
either  experimentally  or  through  another  approximation  method. 

2.  Square  platform  transfer  function.   Consider  plane 
progressive  waves  of  frequency  oj  coming  from  an  angle  3  off  the 
platform  axis.   See  Figure  III-l.   Each  leg  is  thus  subject  to 
waves  of  the  same  frequency  but  different  phase.   For  deep 

water  waves  the  wavenumber  is  k  -  2tt^  ■  a2  .   Thus  the  wave- 

X   g 
length  is 

\  -   22L9L 

The  phase  shift  in  one  wavelength  is  2tt  .   The  phase 
shift,  $ . ,  due  to  a  distance,  1.,  along  the  axis  of  prop- 


agation is  thus 


#.  -  11k  -  si  i± 


Consider  the  center  of  the  platform  as  the  point  of  zero 
phase  shift.   Let  d  be  the  leg  spacing  and  a  =  d//2  be  the 
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distance  from  the  center  of  the  platform  to  the  center  of  any 

leg. 

The  distance  causing  a  phase  shift  for  leg  1  is 

1,  -  a  cos(-  -  3) .   For  leg  2  it  is  1-  =  -a  sin(?  -  3) . 
14  £.  \ 

Similarly  for  legs  3  and  4  the  lengths  are  1_  »  a  sin(-  -  3) 

J         4 

and  1.  -  -a  cos (7  -  3).   From  Equation  III-l  the  phase 
4  4 

shifts  are:  ^2      v 


4,  =  2!  a  cosC^  -  3) 
1  g      4 


6  — S*a  sin(5  -  3) 
a   g       4 


CO2         ,TT     fl» 

<f>  a  -  a  sinv-  -  3) 
3   a       4 


2 

CO  iff  o\ 

a  =--  a  cos(-  -  3) 
-  g       4 


III-2 


Note  however  that: 


cos(-  -  3)  -  cos(?)cos(3)  +  sin(5)sin(3) 
4  4  * 

-  (/2/2) (cos(3)  +  sin(3) } 


sin(?  -  8)  -  sin(?)cos(3)  -  cos(?)sin(3) 
4  4  4 

■  (/2/2){cos(3)  -  sin(3)) 


Letting  T  =  o>2d/g  and  recalling  a  =  (/2/2)d  then 
Equations  III-2  become: 
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T 
<J>  -  -  {cos(P)  +  sin($)  } 
i    2 


4   — -  (cos(3)  -  sin(S)} 

2     2 


=»  ?  (cos(S)  -  sin(3)> 
3    2 


*  »-|  (cos(3)  +  sin(8)>  III-3 

%    2 


3.   Heave.   The  transfer  function  for  heave  for  the 
total  structure,  Tnt,  is  the  sum  of  the  contributions  of 
each  leg.   Taking  phase  into  consideration  each  leg  contri- 
butes Ts  exp{i(o)t+<J> . )  }  .   Thus  the  quantity  of  interest,  the 
square  of  the  magnitude  of  T^,  is 

|Tht|2  =  |rg|2  |ei4,1+  e1*2*  ei(,>3  +  ei<h"|  2 . 

Incorporating  the  identity  e1(p  =  cos(<f>)  +  i  sin(<J>)  , 
summing  real  and  imaginary  parts,  and  taking  the  magnitude 
squared  of  the  exponentials  this  becomes 

|Th  |  *  -  |rg|2  ({cosU^+cosUJ+cos^jJ+costM}2 

+  {sin(<J)i)+sin(<)>2)+sin(<|»3)+sin(<|>H)  }2)        III-4 

Note  however  from  Equation  III-3  that  <K  =  -<frx  and 
$2  =  -$3.   Thus  we  have 
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cosUj)    =  cos($J  sinC^)    =   -sin($„) 

cos(<j>2)    ■  cos(<J>3)  sin(<fr2)    =   -sin(<J>s) 

and  Equation  III-4   becomes 

lruJ*   =    I  Ts  |  2    (2   cos(<h)    +   2   cos(6    )}2 
ht  s  3 

|rht|2  =  4  |r  |2  {cos(«jij)  +  cos(<t>3)}2.        III-5 

By  replacing  $i    and  <J>3  with  their  values  from  Equation 
III-3  and  setting  IL  ■  | r.  | 2/| rg| 2  this  becomes 

n     m   4  (cos(-(cos  3  +sin  3  ))  +  cos{|(cos  3  -sin  3  )))2 

T  T  T  T 

R.  ■  4  icos(-  cos  S)cos(=  sin  3)  -  sin(s  cos  3)sin(£  sin  3) 
n  2  2  2  2 

+  cos(r  cos  3)cos(5  sin  3)  +  sin(«j  cos  3)sin(-  sin  3))2 
2  2  z  2 

T  T 

R.  *  4  {2  cos{-  cos  3)cos(-  sin  S)}2  III-6 

"  2  2 

A  form  often  more  suitable  for  computer  work  is  obtained 
by  applying  the  identity  cos2  A  =  *s  +  h   cos  2A  yielding 

|Tht|2  -  4  (1+  cos(T  cos  3))  (1+  cos(T  sin  3))  |r  |2 

III-7 
which  can  be  implemented  by  recalling  that  T=  o>2d/g. 
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4.   Horizontal  force  assumptions.   The  roll  and  pitch 
moments  experienced  by  a  TLP  will  be  dependent  upon  horizon- 
tal as  well  as  vertical  forces.   In  this  report  it  has  been 
assumed  that  the  total  moment  producing  effect  of  the  hori- 
zontal force  is  small  compared  to  the  vertical  forces.   This 
is  perhaps  not  always  the  case  but  is  assumed  in  this  discus- 
sion because,  while  the  horizontal  wave  forces  may  be  large, 
they  are  distributed  along  the  legs,  both  above  and  below  the 
center  of  rotation  of  the  platform.   The  force  on  the  leg 
near  the  surface  is  large  but  with  a  small  moment  arm.   The 
forces  below  the  center  of  rotation  will  be  smaller  but  will 
extend  further  down  the  leg  and  thus  have  a  large  moment  arm, 
producing  a  moment  similar  in  magnitude  but  opposite  in  sign, 
and  will  tend  to  cancel  the  other.   Thus  the  moment  producing 
effects  of  the  horizontal  forces  are  assumed  to  be  signifi- 
cantly smaller  than  those  of  the  vertical  forces.   When  the 
location  of  the  center  of  rotation  and  the  roll  and  pitch 
natural  periods  are  known,  as  they  would  be  for  any  particu- 
lar TLP  design,  the  horizontal  contribution  to  the  roll  and 
pitch  moment  could  be  calculated.   It  can  be  shown  that  the 
phase  angles  encountered  when  adding  up  the  moments  generated 
by  the  horizontal  forces  on  each  of  the  four  legs  would  be 
the  same  as  calculated  previously  for  heave.   Therefore  the 
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results  computed  for  heave  will  serve  as  a  check  on  the  qual- 
itative behavior  of  the  role  of  horizontal  forces  in  the 
generation  of  roll  and  pitch  exciting  forces. 

5.   Pitch  and  roll.   Pitch  and  roll  for  a  square  plat- 
form are  exactly  symmetric.   Thus  only  roll  will  be  discussed 
here  with  the  understanding  that  pitch  could  be  obtained  by 
shifting  the  coordinates  by  90  degrees.   Taking  the  direction 
of  positive  roll  moment  to  be  as  shown  in  Figure  III-l  the 
total  transfer  function  for  roll,  squared,  would  be 


r  J2  -  |r  |2  |  (d/2)ei*1  +  (d/2)ei(',2+(-d/2)ei*3+(-d/2)ei*M2 
rt       s 


=  lrsl2  I'  le^i+e^-ei^-eiH2 


Let  R  ■  |r  .|2/|r  |2.   Incorporating  the  identity 
e1^  =  cos  $  +  i  sin  <t>    ,  summing  the  real  and  imaginary 
parts,  and  taking  the  magnitude  squared  this  becomes 

R*-  ™  I   {(cos  <f>.+  cos  <fr  -  cos  A  -  cos  $  )2 

XT     4  *  *  3  ' 

+  (sin  4)j+   sin  $2-   sin  <J>3-   sin  <\ )  2 }  III-8 

However,    from  equations   III-3   <f>%    ■   -^    and   <f>a    ■   -<J>2    . 
Thus   cos   ^h   *  cos   4>t    t    cos  4>3    ■  cos  <J>2    t    sin  $„   ■  -sin  4^    , 
and   sin  <J>3=   -sin  4»2    .      Therefore   the  cosine   term  in 
Equation  III-8    reduces    to  zero  and   the   sine   term  becomes 
(2   sin  <!>!    +    2   sin  <J>2)2,    leaving 
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R  =  d   (  sin  $!+  sin  <f>2)  . 

Incorporating  the  values  from  Equations  III-3  this  becomes 

T  T 

R  =  d2  (sin{=(cos  8  +sin  B  ) }+sin{-- (cos  3  -sin  B  )))2 

•i  T  T  T  T 

R_  «  dz  {sin(-cos  3) cos (-sin  8)  +  cos (-cos  B)sin(-sin  B) 
r  2  2  2  2 

T  T  t  T 

+  sin(--cos  B)cos(-sin  B)  +  cos(--cos  B)sin(-sin  B) } 

MM  Mm  & 

R  =  d2  {sin(-cos  B)cos(-sin  B)  +  cos(-cos  8)sin(±sin  B) 
r  2  2  2  2 

mm  T  T  ■» 

-  sin(-cos  B)cos(-sin  B)  +  cos(^cos  8)sin(-sin  B) } 
2  2  »  2 

R_  =  4  d2  cos2 (-cos  B)sin2(?sin  B)  III-9 


By  applying  the  identities  cos2 A  =  h   +  h  cos  2A  and 
sin2A  =  ^  -  h  cos  2A  this  yields  a  form  more  suitable  for 
computer  work. 

|Trt|2  =  d2  {1  +  cos(T  cos  B)}(1  -  cos(T  sin  B) )  |rs|2 

111-10 

As  mentioned  previously  pitch  is  simply  rotated  by  90 
degrees.   Thus  cos  B  becomes  sin  8  and  sin  B  becomes  -cos  B  . 
The  transfer  function  for  pitch  for  a  square  platform  is  thus 

|T   I2  =  d2  {1  -  cos(T  cos  8)H1  +  cos(T  sin  B)  }  |r J2. 
pt  s 

III-ll 
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B.   Wave  spectrum. 

The  wave  spectrum  incident  on  the  platform  varies  with 
both  frequency  and  direction  and  is  denoted  S  (co,3#30). 
So  is  the  center  of  an  arc  of  180   from  which  the  waves 
are  coming. 

For  simplicity  of  calculations  it  is  often  assumed,  as 
it  shall  be  in  this  report,  that  the  total  wave  spectrum 
can  undergo  separation  of  variables.   Thus  we  can  express 
the  spectrum  as 


S  (o),S,3ft)  -  S  (cd)  F(S,3) 


where  the  entire  frequency  dependence  is  contained  in 
Sn(aj)  and  the  entire  angular  dependence  is  in  F(3,30). 
This  section  deals  with  approximations  for  the  angular 
dependence ,  F ( 3 , 0  0 ) . 

1.   Cosine  squared  wave  spectrum.   An  often  used  method 
for  describing  the  angular  spreading  of  waves  is  the  cosine 
squared  wave  spectrum  (see  Figure  III-2) .   It  has  the 
advantage  that  it  vanishes  for  values  of  3  equal  to  plus  and 
minus  it/2.   The  formula  for  the  cosine  squared  spectrum  is 

Sn(co,3,30)  =  C  p(3,30)  SnU)  ;   -u/2  <  3-S0  <  it/2 
where  30  is  the  center  of  the  wave  spectrum  and  C  is  a  constant 
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used   to   normalize    to   unity   the   area  between   -tt/2   and   tt/2 

It  is  derived  as  follows: 

ir/2 
1   -   /  C   cos2(0)    dg 

-TT/2 

1  */2 

-   *    2   /   cos2 (3)    dg    =    2    (7) 

C  0  4 


c-  2 

IT 


Thus   the   equation   for   the  cosine   squared   spectrum  is 

S    (w,S,30)    =  ^  s    (W)    cos2(e-3o)      ;    -5  <    (8-B0)    <  5 
n"irn  2-2 

111-12 
This  spectrum  suffers  from  one  major  disadvantage. 
It  cannot  allow  for  variations  in  the  directionality  of  the 
waves.   This  problem  can  be  eased  somewhat  by  considering 
cosine  to  the  fourth,  to  the  sixth,  and  so  forth  but  this 
results  in  only  step  changes  in  the  directionality.   Methods 
have  been  designed  to  provide  continuous  changes  but  they  are 
somewhat  cumbersome  and  comparison  is  not  apparent. 

2.   Gaussian  wave  spectrum.   The  method  used  in  this  re- 
port to  achieve  continuous  changes  in  directionality  is  to 
consider  the  waves  to  be  spread  about  the  central  angle  8   in 
a  Gaussian  fashion  with  standard  deviation  o.   The  reader  not 
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-IT/2 


Figure  III-2   Cosine  Squared  Spreading 


-tt/2  0  (S-6o) 

Figure  III-3   Gaussian  Spreading 


ir/2 
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familiar  with  Gaussian  (or  normal)  distributions  may  wish  to 
refer  to  any  good  probability  text  for  an  explanation  more 
detailed  than  the  following. 

The  probability  distribution  of  a  Gaussian  process  with 
mean  80  an<*  standard  deviation  a  is  described  by 

p(S,30)  -  — - —  exp{-(3-30)  V(2a2)) 

/2iT  a 

A  spectrum  of  waves  defined  by  a  Gaussian  process  would 
thus  be 

SnU,0,Bo>  =  C  p(3,80)  SnU)  ;  -ir/2<B-30<ir/2    111-13 

where  again  C  is  a  constant  to  normalize  to  unity  the  area 
between  -tt/2  and  tt/2.   It  is  derived  as  follows,  referencing 

30  to  zero: 

tt/2 
1  =   /   C  p($,3o)  d8 

-TT/2 
1      U/2 

£  -   /   p(3,So)  d3 
c  -tt/2 

This  integral  is  evaluated  using  the  properties  of  the 
Gaussian  distribution  and  yields 


C  - 


x 


2  uy-i 


where   *(x)  =~4=  /  exp(-t2/2)  dt 
/2ir  -00 
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is  known  as  the  normal  distribution  function  and  is  tabulated 
in  probability  texts  and  math  tables.   Thus  C  is  a  function 
of  the  standard  deviation  and  shall  hereafter  be  denoted  C(a), 

Unlike  the  cosine  squared  spreading  the  Gaussian  distri- 
bution does  not  vanish  for  3  equal  to  plus  and  minus  tt/2. 
This  is  why  C(a)  varies  with  o.   It  must  account  for  the 
cropping  of  the  "tails"  of  the  distribution  for  |3|  >  tt/2 . 

Some  simplifying  approximations  can  be  made  however. 
As  can  be  seen  from  Figure  III-3  when  a  reaches  a  value  of 
2  or  larger  the  angular  spectrum  approaches  a  constant 
value.   That  value  is  1/tt  (in  order  to  achieve  unity  area)  . 
Thus  Equation  111-13  becomes 

Sn(a>,B,80)  ■  -  sn^)  ''   ~-  <  B-B0  <  \  for  a>2. 

For  a  <  1  the  area  in  the  tails  diminishes  rapidly 
leaving  almost  all  of  the  area  between  -ir/2  and  tt/2.   As 
can  be  seen  in  Figure  III-4  C(a)  rapidly  approaches  unity. 
For  a=1.0  01.13.   For  a=0.9  01.09.   So  for  a  <  1 
Equation  111-13  can  be  approximated  by 

s  («,e,e.>  *  s  u)  {-7=L— }  expC^'fp*2}  . 

For  values  of  a  between  1  and  2  the  entire  exact  formula 
must  be  used 
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C(0) 


0.5        1.0        1.5 
a  (radians) 


2.0 


Figure  III-4   Correction  Factor  for  Loss  of  Gaussian  Tails 
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s  (..■.».>  =  v»>  <jTih^H72h>  ^^f^-1    IIX-14 


-(8-0o)2 

M-y== — i    expt 

_)-l   /2tt  a 
2a 

3.   Comparison  of  cosine  squared  with  Gaussian.   Since 
many  people  are  familiar  with  cosine  squared  spreading  it 
is  useful  to  discuss  for  what  values  of  a  the  Gaussian  would 
approximate  the  cosine  squared.   Although  the  two  functions 
are  quite  distinct  there  are  several  possible  methods  to 
compare  them. 

The  normalized  cosine  squared  distribution  has  a  value 
of  2/tt  -    0.637  for  S  ■  89.  For  Gaussian  spreading  this  would 
be  C(a)  /  (/2ir  a)  at  3  =  So  •   If  a  is  assumed  less  than  1 
then  C(o)  »  1.   Equating  the  two  values  yeilds  a   =  0.627. 

When  8  =  30  +  ir/4  the  cosine  squared  distribution  is 
half  of  its  maximum  value.  This  situation  can  be  matched 
with  a  Gaussian  yielding  still  another  "equivalent"  a  for 
the  cosine  squared. 


C(ai)     f-(7r/4)2,       ex  C(m) 


exp{"(!/})  }  -  (0.5) 
2a  z 


°   =  {   ln(0.5)   } 


o  =  0.667 
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The  two  values  obtained  are  approximately  equal.   Thus  the 
Gaussian  distribution  with  a=0.65  is  very  similar  to  the  more 
familiar  cosine  squared  distribution.   A  smaller  a  means  a  more 
unidirectional  wave  front  and  a  larger  a   means  more  spreading. 

C.   Platform  exciting  forces. 

Each  differential  element  of  the  spectrum  S  (co,B,B0)  d8 
induces  mean  square  force  or  moment  upon  the  platform  through 
the  appropriate  platform  modal  force  transfer  function, 
magnitude  squared,  |rt(o),B)|2.   Thus  the  total  force  or 
moment  spectrum  is  given  by 

TT 

H(<o,Bo)  =  /  |T  (&>,B)|2  S  (a),8,30)  dB      111-15 

t         n      " 

—  TT 

It  is  important  to  note  that  this  modal  force  spectrum 
is  dependent  upon  80 ,  the  center  of  the  directional  spectrum. 

1 .   Square  platform,  heave  force.   The  heave  force 
spectrum  for  a  square  platform  in  normally  distributed 
waves  is  derived  by  substituting  Equations  III-7  and  111-14 

into  Equation  111-15  yielding: 

B0+ff/2 
H  (o>,80)  -  /    |T  |2  4{l+cos(T  cos  B)}{l+cos(T  sin  B)  }x 
8o-*/2   S 

S  M    { 1 H-r=L— )  exp{~(^~go2)2}  dB 

n""   2*(U-1  /2¥o    *        2   az 

2a 
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H.U.B.)   =  4   |r J2  s   <«•>   ^f^* 


So+*/2  ,R    fi    ,2 

/        {l+cos(T   cos   6)}(l+cos(T   sin   8)  >    expC"*p~pK    }    d8   111-16 


80-*/2 


2  a^ 


Setting  In(<o,$  )  equal  to  the  integral  in  Equation 
111-16  and  expanding  the  integrand  yields 


30+tt/2 
I    (a),80)    ■        /  (I   +  cos(T  cos   8)    +   cos(T   sin   3)    + 

h  3, -ir/2 

cos(T  cos  3)cos(T  sin  8))  exp{"(^ffl)  }  d8    111-17 


where  H.  <«,B0)    =   4    |r    |2    S_(g>)    £J£1     I    (»,30>    . 
n  °  s'         n  /2tt   a      h 

A  closed  form  solution  of  Equation  111-17  does  not 
exist.   Thus  a  numerical  integration  was  necessary 
(  see  Section  IV  ) . 

2.   Square  platform,  roll  moment.   Similarly  the  roll 
moment  spectrum  for  a  square  platform  in  normally  distributed 
waves  is  derived  by  substituting  Equations  111-10  and  111-14 
into  Equation  111-15  yielding: 
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B0+t/2 
H    (o),60)    =   /  d2    (l+cos(T   cos    S)}(l-cos(T    sin   3)  }    |r    |  2* 

r  S.-1T/2  S 


3(a))     {      J    ■       )    {7=L-}    exp{-(S-8o)2} 
n  2    $  (I_J-1        /2tt   a  2 

2a 


2 

dB 


0)  =  d2   |r  |2  s  (co)  £&- 

1    s'         n  /2ir   a 


B„W2  ,g    s    v  2 

/         {l+cos(T   cos    3)}{l-cos(T   sin   3)}    exp{    ~    2         }    d3    111-18 

So-TT/2  *° 


Setting  I  (<u,30)  equal  to  the  integral  in  Equation  111-18 

and  expanding  the  integrand  yields 

S0+tt/2 
IrU,0o)  ■   /   {l+cos(T  cos  3)-cos(T  sin  3)- 
3o-tt/2 


cos(T  cos  3)cos(T  sin  3)}  exp{~(jj~?o)  )  d3     111-19 

2c 


where  Hr<a»,B.)  =  d2  |TS|2  Sn(o»)  ^|JJ£L  Ir(«,Bo)  . 


Again  a  closed  form  solution  does  not  exist  and  the  same 
numerical  integration  technique  was  applied  (  see  Section  IV) . 
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3.   Square  platform,  pitch  moment.   As  mentioned 
previously,  for  a  square  TLP  the  pitch  moment  is  simply 
rotated  90  degrees  from  the  roll  moment.   Thus 


h  <«,&,)  =  d2   |rs|2  s  («)     cjoL  i  <«,$0) 

v  /2ir   a      v 


where      I    (a>,B0)    ■  I    (a>,B0+T/2)     . 

ST  * 

4.   Normalization.   In  order  to  compare  the  effect  of 

the  wave  spreading  for  different  configurations  a  method 

of  normalization  was  developed.   The  normalization  factor 

was  taken  to  be  the  force  or  moment  frequency  spectrum 

resulting  from  waves  of  equal  amplitude  from  all  directions 

in  one  half  of  the  plane.   Thus  Sn(o>,B,6o)  equals  (1/tt)  Sn(oj) 

for  (Bo-ff/2)  <  B  <  (Bo+t/2)  and  zero  elsewhere.   The 

normalization  frequency  spectrum  becomes 

$o+^/2 
H    U,Bo)    -      /  |T    <o),B)|2     (1/tt)    S    U)    dB  . 

°  Bo-tt/2 

Because  of  the  symmetry  of  a  square  platform  and 

the  uniformity  of  this  spreading  it  makes  no  difference 

where  the  center  of  the  spectrum  is.   For  convenience 

B0  was  taken  to  be  zero  and  thus  the  normalization  spectrum 

becomes 

1       w/2 
H0(co)  =  Z   Sn((0)  f      |r.(Tr,B)|2  dB        IH-20 

-tt/2 

30 


which  is  simply  the  mean  square  value  of  T      times  S  (a))  . 

t        n 

For  heave  and  roll  they  are  as  follows 


TT/2 


HOh((u)  =  ff   Sn(<u)     'rs'         f    *1+cos<T  cos   3)Hl+cos(T   sin   3))    d3 

111-21 

,2  w/2 

Hrt  ((u)=2  s  (to)  |r  |2   /  (l+cos(T  cos  3)Hl-cos(T  sin  3))  d3 
°r    *   n      S   -*/2 

111-22 

The  normalized  spectrum  for  heave  is  thus 

H,  U,3o) 

HN  (co,30)  =  -2 — -  *  111-23 

h     °      H-.  (<o) 
Oh 

Expanding  the  numerator  and  denominator  and  cancelling 

the  common  terms  S  (<o)  and  |r  I2  this  leaves 

n  s 

m    (<,,S0)  =  1^<"'f°>  111-24 

h  I3n(w) 

where 

ti  /   q  ^     Hh(<a,gB) 

I1h(tt'8«)  =  S  (o.)  |r  |2 
n     '  s' 

4   C(a)   e°+7r/2 
3  7==  ^      /    (l+cos(T  cos  3)}  * 

/I?  a    B,-ir/2 


(l+cos(T  sin  3)}  exp{~ (B~?o } * }  d3 

2a2 


13    (to)    -  -Ok 


h'  S    (a>)    |r 

n  '    i 


4        ff/2 

■    (-)      /    {l+cos(T  cos   3)Hl+cos(T   sin   3))    d3 
-ir/2 
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Similarly  for  roll 


Hr(co,30)    II  (u,B0) 

HNr(^S»)  =   OB   =   13  M 
Or  r 


where 

Hr(o>,3o) 


II  (o>,30)  - 


r        Sn(o,)  |rs|2 


3o+tt/2 
=  _A_  c<q>    /    (l+cos(T  COS  S)}  X 
/2?  a   30-tt/2 


(l-cos(T   sin   3)}    expT^'f0*    }    d3 


I3rU)  '  vSfTr^ 


d2  */2 


/    {l+cos(T  cos   3)}{l-cos(T   sin   3)  )    d3    . 
"  -ir/2 
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IV.  Numerical  Prediction  of  Exciting  Forces 

A.  Discussion. 

Because  of  the  complicated  forms  of  the  integrands  of 
Equations  111-17  and  111-19  no  closed  form  solutions  to 
the  integrals  were  found  to  exist.   Numerical  methods  were 
employed  to  determine  the  variation  of  resulting  forces 
caused  by  the  variation  in  input  parameters.   Appendices 
I  and  II  provide  program  listings  and  an  explanation  for 
their  use. 

B.  Assumed  values. 

In  order  to  use  numerical  methods  of  calculation 
values  must  be  assumed  for  the  relevant  parameters.   In 
this  thesis  the  initial  values  assumed  come  from  two 
Amoco  papers  published  in  the  "Journal  of  Petroleum 
Technology."   (See  Reference  1.)   These  papers  discuss 
a  proposed  vertically  moored  platform  (VMP)  which  is 
identical  to  the  concept  of  a  TLP.   Amoco' s  VMP  is  proposed 
to  have  leg  spacings  of  160  feet  and  is  estimated  to  have 
heave,  pitch,  and  roll  natural  periods  of  approximately  two 
seconds.   Because  of  the  popularity  of  the  cosine  squared 
wave  spectrum  a  comparable  standard  deviation  of  0.6  5 
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radians  was  chosen. 

The  parameter  of  interest  in  all  of  the  following 
analyses  is  HN,  the  ratio  of  the  magnitude  of  the  transfer 
function  for  Gaussian  wave  spreading  to  the  magnitude  of 
the  transfer  function  for  uniform  wave  spreading. 

When  the  Amoco  TLP  was  analyzed  in  normally  distributed 
waves  with  a  0.65  standard  deviation  the  results  were  some- 
what surprising.   The  variation  in  HN  from  unity  in  heave 
was  less  than  one  percent  for  all  values  of  wave  spectrum 
incidence,  $0 ,    (see  Table  IV-1) .   A  much  greater  variation 
of  HN  for  roll  was  expected  but  HN  still  remained  within 
ten  percent  of  unity  (see  Table  IV-2) .   Clearly  for  this 
case  an  initial  assumption  of  uniform  wave  spreading  and 
the  resulting  simplifications  in  calculation  of  the  force 
spectrum  will  provide  very  useful  first  pass  information 
to  the  designer. 

However  the  Amoco  proposal  is  by  no  means  the  only 
viable  configuration  for  a  TLP.   One  must  consider  the 
variation  in  HN  as  other  parameters  vary. 

c«   Variation  in  leg  spacing. 

The  most  obvious  variation  in  a  platform  is  the  spacing 
between  the  legs.   With  the  160  feet  of  the  Amoco  platform 
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as  a  reference,  leg  spacings  ranging  from  100  to  250  feet 
were  considered.   Analysis  results  for  heave  are  presented 
in  Table  IV-3  and  for  roll  (and  hence  pitch)  in  Table  IV-4. 
For  all  of  these  configurations  HN  for  heave  is  still 
within  one  percent  of  unity.   Roll  suffers  somewhat  but 
the  calculated  HN  stays  well  within  20%  of  one.   For  initial 
design  work  this  accuracy  is  sufficient  and  the  assumption 
of  uniform  spreading  would  be  justifiable. 

D.   Variation  in  period. 

Because  the  wave  spectrum  is  continuous  in  frequency 
waves  of  all  periods  will  be  incident  upon  the  structure. 
However,  the  period  of  paramount  concern  is  the  natural 
period  of  the  structure.   This  varies  from  structure  to 
structure  and  is  somewhat  controllable.   Again  taking  the 
Amoco  value,  two  seconds,  as  a  reference,  periods  from  one 
to  three  seconds  were  considered.   The  leg  spacing  was 
assumed  to  be  160  feet  and  the  standard  deviation  0.6  5 
radians.   Analysis  results  for  heave  are  tabulated  in 
Table  IV-5.   The  results  for  roll  however  vary  enough  to 
be  visible  on  a  polar  plot  of  HN  vs.  30 /  the  center  of 
the  wave  spectrum,  and  are  shown  in  Figure  IV-6.   When 
possible  a  plot  has  been  preferred  to  a  table  because  of 
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the  greater  amount  of  information  compactly  presented. 

For  heave  the  variation  in  HN  remains  at  about  one  per- 
cent reaching  a  high  of  two  percent  for  a  three  second  per- 
iod.  For  roll  HN  shows  marked  fluctuation  but  again  stays 
within  twenty  percent  of  unity.   While  the  natural  period 
has  significant  implications  with  respect  to  the  periods  of 
significant  wave  energy  it  does  not  appear  to  significantly 
alter  the  modal  force  transfer  functions  as  long  as  some 
directional  spreading  is  present. 

E.   Variation  in  standard  diviation 

The  final  parameter  to  vary  is  the  standard  deviation. 
The  standard  deviation  for  uniform  wave  spreading  is  infinite 
while  that  for  unidirectional  waves  is  zero.   However,  as 
shown  in  Section  III-B  standard  deviations  above  two  radians 
approximate  the  condition  of  uniform  wave  spreading.   There- 
fore almost  all  of  the  effect  of  varying  standard  deviation 
occurs  between  zero  and  two.   One  would  thus  expect  a  signi- 
ficant change  in  HN  as  standard  deviation  changes,  particu- 
larly for  small  values. 

This  expectation  is  correct.   Standard  deviation  was 
varied  from  0.125  to  2  radians  for  both  heave  and  roll 
(D  ■  160  feet,  Per  ■  2.0  seconds) .   The  results  are  plotted 
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in  Figures  IV- 7  and  IV-8  for  heave  and  roll  respectively. 

Consider  heave  (Figure  IV-7) .   For  standard  deviations 
of  1.0  or  2.0  radians  HN  stays  within  one  percent  of  unity. 
When  a  =  0.5  the  variation  becomes  five  percent.   For 
a  ■  0.25  it  exceeds  twenty  percent  and  when  the  standard 
deviation  is  0.125  radians  even  the  pattern  of  the  curve 
changes.   For  roll  (Figure  IV-8)  the  results  are  similar, 
although  somewhat  more  pronounced. 

However,  if  one  can  allow  a  15%  error  in  an  initial  de- 
sign approximation  then  the  uniform  wave  spreading  assumption 
is  valid  for  wave  directional  spectra  with  standard  deviations 
greater  than  approximately  0.5  radian.   Note  that  cosine 
squared  spreading  is  included  in  this  category. 

F.   Summary 

Numerous  other  cases  were  considered,  simultaneously 
varying  leg  spacing,  period,  and  standard  deviation.   In  all 
cases  where  a  >  0.5  the  resulting  normalized  amplitude- 
platform  force  transfer  function  had  a  magnitude  within 
20%  of  unity.   Thus  for  waves  with  any  reasonable  amount  of 
spreading  sufficient  accuracy  for  initial  design  calculations 
can  be  obtained  by  using  the  mean  square  value  of  r(a>,6). 
This  result  allows  for  a  simpler  initial  calculation  of 
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platform  response  to  random  waves  in  the  manner  of 

4 
Vandiver.    For  situations  when  a    <  0.5  or  when  more 

accuracy  is  required  Vandiver' s  method  can  still  be 

employed  using  the  computer  programs  of  this  thesis  to 

determine  HN  (Vandiver 's  CI). 
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SQUARE  TENSION  LEG  PLATFORM,  HEAVE  EXCITATION 


WITH  GAUSSIAN  WAVE  DISTRIBUTION 


LEG  SPACING  - 

PERIOD  - 

STANDARD  DEVIATION  = 


160   feet 
2.0   seconds 
0.6  5  radians 


Correction  factor  due  to  loss  of  tails  ■  1.016 
Factor  due  to  uniform  spreading  (13)  ■  3.9  4 

e0  ii 

0  3.90 

5  3.91 

10  3.92 

15  3.93 

20  3.94 

25  3.95 

30  3.96 

35  3.96 

40  3.96 

45  3.97 


HNh  = 

=11/13 

0, 

.990 

0. 

,992 

0, 

.994 

0, 

.997 

1, 

.000 

1, 

.002 

1, 

.004 

1, 

.006 

1. 

.006 

1, 

,007 

Table  IV-1  Normalized  Platform  Heave  Exciting  Force 
(D«160  feet,  Per=2.0  sec,  a=0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  ROLL  EXCITATION 


WITH  GAUSSIAN  WAVE  DISTRIBUTION 


80 

0 
5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 


LEG  SPACING  -        160   feet 
PERIOD  ■  2.0   seconds 

STANDARD  DEVIATION  =0.65  radians 

Correction  factor  due  to  loss  of  tails  ■  1.016 
Factor  due  to  uniform  spreading  (13)  =  0.235  E+05 

II 


0.216 
0.216 
0.217 
0.218 
0.219 
0.222 
0.224 
0.227 
0.231 
0.234 
0.238 
0.241 
0.244 
0.248 
0.251 
0.254 
0.255 
0.256 
0.257 


E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+0  5 
E+05 


H»r» 

=11/13 

0, 

.919 

0, 

,920 

0. 

.922 

0, 

.925 

0. 

.933 

0, 

.943 

0. 

.954 

0, 

.967 

0. 

.981 

0. 

.996 

1, 

.011 

1. 

.026 

1. 

.040 

1, 

.053 

1, 

.068 

1. 

.079 

1. 

.086 

1, 

.090 

1. 

.092 

Table  IV-2  Normalized  Platform  Roll  Exciting  Moment 
(D=160  feet,  Per=2.0  sec,  a=0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  HEAVE  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 


PERIOD  ■  2.0   seconds 

STANDARD  DEVIATION  ■   0.65  radians 

CORRECTION  FACTOR  DUE  TO  LOSS  OF  TAILS  =  1.016 

HNh 
(max.) 

1.003 
1.003 
1.006 
1.004 


LEG  SPACING 

HNh 

(feet) 

(min. 

100 

0.997 

150 

0.998 

200 

0.994 

250 

0.995 

Table  IV-3  Normalized  Platform  Heave  Exciting  Force 
(Various  leg  spacings,  Per=2.0  sec,  0=0.6  5  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  ROLL  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 

PERIOD  =  2.0   seconds 

STANDARD  DEVIATION  =   0.65  radians 

CORRECTION  FACTOR  DUE  TO  LOSS  OF  TAILS  ■  1.016 

HNr 
(max.) 

1.027 
1.067 
1.114 
1.160 


LEG  SPACING 

HNr 

(feet) 

(min.) 

100 

0.978 

150 

0.933 

200 

0.882 

250 

0.846 

Table  IV-4  Normalized  Platform  Roll  Exciting  Moment 
(Various  leg  spacings,  Per=2.0  sec,  a«0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  HEAVE  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 


LEG  SPACING  ■         160   feet 
STANDARD  DEVIATION  ■   0.65  radians 

CORRECTION  FACTOR  DUE  TO  LOSS  OF  TAILS  -  1.016 


PERIOD 

(seconds) 

1. 

,0 

1. 

,5 

2. 

,0 

2. 

,5 

3, 

,0 

HN.  HNV 

(mm.)  (max.) 


0.997  1.001 

0.996  1.005 

0.990  1.007 

0.994  1.004 

0.980  1.024 


Table  IV-5  Normalized  Platform  Heave  Exciting  Force 
(D=160,  Various  periods,  a=0.65  rad.) 
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Figure  IV-6   Normalized  Platform  Roll  Exciting  Moment 
(D=160,  various  periods,  c=0.6  5  rad.) 
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B0  =  90' 


-  0.125 


HN^  ■ 


Hh(ojfe0) 

H0h("> 


2.0  M° 


Bfl  =  0 


Figure  IV-7   Normalized  Platform  Heave  Exciting  Force 
(D=160,  Per  =2.0  sec,  various  a) 
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B0  =  90 


HN   = 


Hr(a),g0) 
H0r(co) 


=0' 


Figure  IV-8  Normalized  Platform  Roll  Exciting  Moment 
(D«160,  Per  ■  2.0  sec,  various  a) 
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V.   Suggestions  for  Experimental  Verification 

A.  Single  leg. 

Because  the  results  of  the  computer  program  provide  a 
transfer  function  between  forces  on  a  single  leg  to  platform 
forces  and  moments,  the  single  leg  forces  must  be  obtained. 
This  can  be  done  theoretically  or  experimentally.   This 
report  will  deal  with  experimental  determination. 

The  single  leg  is  assumed  to  be  of  circular  cross  sec- 
tion.  Thus  there  is  no  angular  dependence  and  testing  can 
be  performed  in  a  ship  model  towing  tank  with  a  wave  gener- 
ator.  In  the  proposed  test  the  single  leg  model  will  be  held 
rigid  with  load  cells  measuring  the  forces  as  a  result  of 
incident  waves.   These  forces  will  be  normalized  to  forces 
due  to  scaled  unit  amplitude  waves  using  a  linear  assumption 
for  the  force  to  wave  amplitude  relation.   By  performing 
this  test  with  a  variety  of  wave  frequencies  a  plot  of  rs(oj) 
can  be  obtained.   From  this  the  computer  will  predict 
platform  heave  force  and  roll  moment  transfer  functions. 

B.  Square  TLP. 

To  completely  test  the  results  of  this  study  would  re- 
quire  the  ability  to  generate  a  directional  wave  spectrum. 
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This  ability  is  not  available.   However,  the  total  transfer 
function,  r  (a),3),  can  be  experimentally  verified  for  heave, 
pitch,  and  roll. 

Once  the  results  for  a  single  leg  are  obtained  a 
similar  experiment  with  the  total  platform  can  be  performed. 
The  model  should  be  held  rigid  with  load  cells  measuring 
forces  in  all  six  degrees  of  freedom.   The  direction  of 
wave  incidence  can  be  varied  by  turning  the  model. 

This  experiment  will  generate  a  family  of  curves  of 
r.(a),3)  ,  one  curve  for  heave  and  one  for  roll  for  each  angle 
of  incidence  3.   These  curves  can  then  be  compared  to  those 
predicted  using  the  measured  forces  for  rs(w)  and  the  heave 
and  roll  algorithms  developed  in  this  paper,  Equations  III-7 
and  111-10. 
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VI.   Conclusions 

A.   Assumptions  and  approximations. 

1.  Linearity.   The  forces  computed  in  this  paper  are 
those  exerted  by  waves  on  a  motionless  structure.   The 
forces  are  assumed  to  be  linearly  related  to  wave  amplitude. 
Furthermore  the  forces  exerted  on  a  motionless  structure 
can  only  be  useful  in  predicting  dynamic  response  if  the 
equations  of  motion  are  linearized.   For  a  TLP  the  equations 
of  motion  will  be  linear  for  low  to  moderate  sea  states. 
This  is  due  to  the  large  dimensions  of  the  major  members  and 
to  the  relatively  small  response  amplitudes.   Damping  terms 
for  the  TLP  may  also  be  linearized  because  the  damping 
ratios  are  expected  to  be  small  in  heave,  pitch,  and  roll. 

2.  Effects  of  cross  members.   The  model  of  a  TLP  used 
to  derive  the  transfer  functions  had  four  legs  acted  upon  by 
wave  forces,  and  nothing  else.   Real  TLP ' s  must  have  cross 
members  between  the  floatation  legs  to  provide  structural 
support.   These  would  also  be  subject  to  wave  forces. 

The  diameter  of  a  cross  member  will  be  much  smaller  than 
that  of  the  leg.  In  the  Amoco  platform  referred  to  in  Refer- 
ence 1  the  cross  members  have  an  outer  diameter  of  5.5  feet 
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while  the  floatation  leg's  diameter  is  30  feet  near  the  sur- 
face where  the  large  wave  forces  are.   This  is  a  factor  of 
about  5*5  difference.   The  displaced  volume  is  proportional 
to  the  cross  sectional  areas  which  are  almost  a  factor  of  30 
different.   In  addition  the  support  members  are  mostly  near 
the  bottom  of  the  legs,  away  from  the  largest  wave  forces. 
The  parts  that  are  shallow  are  close  to  the  legs  and  would 
react  in  phase  with  them.   Thus  the  effect  of  the  waves  on 
the  cross  members  can  be  assumed  to  be  negligible  with 
respect  to  the  effect  on  the  legs. 

3 .   Wave  disturbances  of  the  legs.   Another  aspect  not 
considered  in  this  model  is  the  effect  of  one  leg  upon  the 
waves  incident  upon  the  others.   In  the  Amoco  example  the 
outer  edges  of  the  legs  at  the  water  surface  are  130  feet 
apart  which  is  only  about  four  times  their  30  foot  diameters 
There  could  be  some  shielding  or  reflecting  effects  present, 
causing  the  results  of  this  study  to  be  in  error.   As  the 
ratio  of  leg  spacing  to  leg  diameter  increases  these  errors 
will  decrease. 

The  primary  purpose  of  this  thesis  was  to  evaluate  the 
influence  of  wave  spreading  on  the  heave,  pitch,  and  roll 
force  spectra.   The  most  important  geometric  effect  is  the 
relative  phases  of  forces  exerted  on  the  legs.   In  cases 
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that  the  wave  spreading  is  shown  to  be  sufficient  to  average 
out  the  relative  phase  effects,  then  it  can  certainly  be 
argued  that  interaction  effects  between  legs  would  also 
be  suppressed. 

B.   Utility  of  results. 

1.  Initial  design  approximation.   For  wave  directional 
spreading  with  a  standard  deviation  of  greater  than  0.5 
radians  the  heave  force  or  roll  or  pitch  moment  imposed 
upon  a  TLP  can  be  approximated  within  20%  to  be  that  imposed 

by  uniformly  spread  seas.   By  incorporating  this  into  the 

4 

method  used  by  Vandiver  an  initial  estimate  of  the  response 

of  a  TLP  can  be  obtained  which  would  be  applicable  in 
preliminary  design. 

2.  More  accurate  approximation.   If  the  standard 
deviation  of  the  wave  spectrum  is  less  than  0.5  radian  or 
if  an  answer  with  more  than  20%  accuracy  is  required  then 
the  computer  program  outlined  in  Appendices  I  and  II  can 
be  used  to  determine  HN  (Vandiver1 s  CI).   This  can  then 

be  incorporated  into  the  calculations  of  the  modal  response 
of  the  structure. 

3.  Variation  to  fit  other  structures.   Results  for 
other  square  TLP ' s  can  easily  be  obtained  simply  by  changing 
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the  values  assigned  for  the  leg  spacing  and  the  period  of 
interest.   If  the  structure  of  interest  is  not  square  the 
same  concept  can  be  applied  by  changing  the  program  defini- 
tions of  the  transfer  functions.   In  this  manner  results 
can  be  obtained  for  any  shape  of  TLP  (square,  triangular, 
rectangular,  etc.)  with  any  leg  spacing  and  any  natural 
period. 

4.   Wave  directional  spectrum  variation.   The  greatest 
advantage  of  modeling  the  wave  directional  spectrum  as  a 
Gaussian  process  is  the  continuous  variation  in  spreading 
it  provides.   The  standard  deviation  of  the  wave  spectrum 
at  a  proposed  location  for  a  TLP  could  be  measured  under 
various  conditions,  noting  that  the  periods  of  the  signi- 
ficant wave  energy  would  also  change.   This  information, 
standard  deviation  and  period,  could  then  be  input  to  this 
program  to  estimate  platform  forces  and  moments  under  those 
conditions.   The  final  information  of  forces,  moments,  and 
responses  would  then  be  applied  to  sizing  the  tension 
members  for  adequate  tensile  and  fatigue  strength  and  to 
ensure  the  platform  motions  do  not  cause  equipment  malfunc- 
tions or  personnel  discomfort. 
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NOMENCLATURE 

a    distance  from  center  of  TLP  to  center  of  any  leg 

C    constant  to  normalize  the  area  of  the  wave 
directional  spectrum 

d    leg  spacing  of  TLP 

g    acceleration  of  gravity 

h    as  a  subscript,  denotes  heave 

H    frequency  spectrum  of  the  total  force  or  moment 
on  the  TLP  due  to  wave  spreading 

Hn    frequency  spectrum  of  the  total  force  or  moment 

on  the  TLP  due  to  uniformly  distributed  waves 

HN    normalized  frequency  spectrum  of  the  total  force 
or  moment  on  the  TLP 

I     integral  portion  of  H   dependent  on  wave  spreading 

Ilx   numerator  of  HN  ,  dependent  on  platform  geometry 
and  wave  spreading 

13    denominator  of  HN  ,  dependent  only  on  platform 
geometry 

k    wavenumber 

1.    distance  from  center  of  TLP  to  projection  of  leg  i 
on  the  direction  of  wave  propagation 

p    as  a  subscript,  denotes  pitch 

r    as  a  subscript,  denotes  roll 

R    ratio  of  magnitude  squared  of  platform  T  to 
magnitude  squared  of  single  leg  r 

s    as  a  subscript,  denotes  single  leg 

Sn   wave  spectrum 
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t    as  a  subscript,  denotes  total  platform 

T    non-dimensional  quantity  determined  by  leg  spacing 
and  wave  frequency 

(3    angle  of  wave  incidence  measured  from  roll  axis 

80    angle  of  incidence  of  center  of  wave  spectrum 

T    transfer  function  from  incident  wave  amplitude  to 
x 

force  (or  moment)  exerted  upon  the  structure 

6-^  Dirac  delta  function 
mn 

A    wavelength  of  incident  waves 

a  standard  deviation  of  wave  directional  spectrum 

£    summation 

<f>.    phase  shift  of  the  wave  incident  on  leg  i  as  compared 
to  the  wave  at  the  center  of  the  TLP 

$    normal  distribution  function 

a)    radian  frequency  of  incident  waves 
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Appendix  I .   Computer  Program  Listing 

A.   Discussion. 

Because  of  the  complicated  forms  of  the  integrands  of 
Equations  111-17  and  111-19  no  closed  form  solution  to  the 
integral  was  found  to  exist.   Yet  the  trigonometric  func- 
tions with  trigonometric  arguments  caused  such  rapid  oscil- 
lations of  the  integrand  that  most  numerical  methods  failed 
to  converge.   The  method  found  to  be  the  most  useful  was 
Gauss-Legendre  quadrature  (see  Appendix  III) .   A  quadrature 
of  fifth  order  was  found  to  give  sufficient  numerical 
accuracy  and  is  contained  in  the  subroutine  GAULEG. 

Because  the  wave  directional  spectrum  was  only  non- 
zero between  plus  and  minus  it/2  radians  a  correction  factor 
was  necessary  to  ensure  its  area  was  unity.   This  correction 
was  accomplished  by  an  eighth  order  polynomial  approximation 
to  the  area  under  the  normal  curve.   This  approximation  was 
developed  by  the  author  and  is  contained  in  the  subroutine 
ERF. 

Because  of  peculiarities  in  the  computer  system 
available  the  platform  transfer  functions  are  entered  as 
function  subprograms  HEAV  and  FHEAVE  for  heave,  and  ROLL 
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and  FROLL  for  roll,  the  former  of  each  pair  being  the 
transfer  function  with  uniform  wave  spreading. 

The  program  is  written  in  FORTRAN-IV  language  and  has 
numerous  comment  statements  to  aid  in  its  understanding. 

B.   Program  listing. 

A  listing  of  the  programs  for  heave  and  roll,  with 
associated  subprograms,  follows. 
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Z  TFSSTDS  LF3  PLATFOR* 

C  HF1  VF 

EXTERNAL  WEAV^Ff'SHV 

C  SET  STANDARD  DEVIATION"  (SD),LF3  SPACES  (P),AVD  PERIOD  (PER) 
D»152. 

S0=0.55 
PrR=2. 

z   calculate  v»riakce  a*»0  t 
200    y\?=;d**2 

300     T=1.?2SC3R*r>/»rR**2 
«T-  (5,5) 

5    foptat  (m',n?,  'square  tf^isiov  le3  pljtfor", reave  excitatids') 
writ:  (5,1?) 

10    FORMAT  C  ,rri5,,3A033IA'«  WAVE  DISTRI BOTI nv//) 

WRITE  (5,15)  D 
15    FORMAT  (•  'jTIS^LEG  5?  5  cn'~='  ,  "UD,  F5.0  ,  '  FEET') 

HPITE  (5,20)  Prt> 
20     FORMAT  C  *,rie,  'PEPIODs'  ,TUO,**.1  ,  '  SSFONO^') 

WRITE  (5,25)  3D 
25     FOR*AI  C  ',T1C,  'STANDARD  DEVIATIOVs*  ,T33  ,Ff  .2,'  RADIANS'//) 

-  COPRECTION"  FACna  FDR  LOSS  OF  TAILr  OF  SAUSSIAR 

*TOTH=1  .*7:i795/SD 

CALL  FP.F  (*I0rH,P-lI) 

C=1./(2.*PUI-1.) 

?'!"  (5,33)  : 
30  FOF'"^  C  ',T15,' CORRECTION  FACTOR  OOF  ro  LOSS  OF  TAILS=* ,F6.3) 
3  FACTOR  OUe  TO  UNIFORM  SDR??.DT  "^ 

YL=-1. 570756 
ICO*  1.570796 
D  P  =  0 . 

CALI  "ArTT."T  (<L,X'J,HrAV,R,I",n'»,VAP,T) 
RI3=R 

WIT"    (5,35)    RI3 
35  FO?w«T    C     *,T15, 'FACIOR    DOE    'D    TJ1IFOPW    FSRF«DT1S    (I3)  =  ',E11.U) 

KRITF    (5,!»D    IP 

uo        fop.'mt  r    •,r2r, '(frro0  covririo's',13,')  •/) 

w'RITF.    (5,U5) 
u5    FOP*A?    ('     ','    BET*.'  ,T23,  '  TV  ,  'us,  '  KS  =  I1/I3 ' , T52,  'ERROR?'/) 
C 

Z    VA^Y    DIDrCrT0'J    OF    CENTER    Or    iP.VF    i^^rTPtJ* 
PODESsO. 
100.         SO=B0DE5*1.7'45323E-2 
:<L1=KL+9<? 

CALL    3AULE3    ( ?L ", XU 1 ,FRFJ 7 ,Y ,IEP,PO ,VAP,T) 
?71^Y*C/F0 
4:T  =  RT1/PI3 

WRITE    (6,50)    B0*F.1,RT1,«"J,TF.P 
50    FOR*AT    (.     ',F5.1  ,-l7,S11.ttfTU2,FlO.,4,T52,n) 
P0r>-l  =  B^Dr-+5. 

IF    (BODES    .LF.    "5.)    f*3    TO    100 
STOP 
E?ID 
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C  TENSION  LEG  PLATFORM 

C  ROLL 

EXTERNAL  ROLL.FRCLL 
C  SET  STAMOARO  DEVIATION  (SD).LEG  SPACING  (G)»ANO  PERIOD  (PER) 
0=160. 
SD  =  C65 
PER=2. 
C  CALCULATE  VARIANCE  AND  T 
2(C      VAP=SD**2 
300     Trl.22603fi*C/PER«*2 
WRITE  (6,5) 
5     FORMAT  (  *l*»T10.*Sr,U.'RE  TENSION  LEG  PLATFORM.  =(OLL  EXCITATION*) 
WRITE  (6.10) 
10     FORMAT  (»  *.T15. 'GAUSSIAN  WAVE  DISTRIBUTION*//) 

WRITE  (6. 1  "5)  0 
15    FOovAT  <♦  *.T15.»LEG  SF  ACINGs  •  ,T<»0.F5.  0  t  •  FEET") 

WRITE  (6.2H)  PER 
2C     FORMAT  (•  •♦T15«»PERIOE=». Tt0.F5.lt*  SECONDS*) 

WRITE  (6.25)  SD 
25     FORMAT  (*  * »T15* 'STANDARD  DE V  IA T  I  ON= * t T39 ,F6. 2 . *  RADIANS*//) 
C 

C  CORRECTION  FACTOR  FOR  LOSS  OF  TAILS  CF  GAUSSIAN 

WIDTH=1.570796/S0 
CALL  ERF  (WIDTH.PHI) 
C=1./C2.*PHI-1.) 
WRITE  (6.30)  C 
30  FORMAT  (*  • , T15» •CORRECTI ON  FACTOR  DUE  TO  LOSS  OF  T A  I LS  =  * *  Ft. 3 ) 
C  FACTOR  DUE  TO  UNIFORM  SPREADING 
C 

XL=-1. 570796 
XU=  1.570796 
E"  =  0. 

call  qaulen  (xl*xu»roll.r«ir.r-0»VAP»T) 
RI3=R*D»*2 
WRITE  (6.35)  PI3 
35     FORMAT  (•  *  .T15 . *FAC TOR  DUE  TO  UNIFORM  SPREADING  < I  3 )  =  » . El 1 . 4 ) 

WRITE  (6.40)  1R 
♦  0      FORMAT  (•  •  »T20»* (ERPOR  COND I T I ON= * . I  3 /) 

WRITE  (6.45) 
4*  FORMAT  (*  ».*  PETA  *.T27  ♦•  I1*»T45  »  *HNrH/I  7,*t  T62.  *:SR0R?*/) 


C  VA9Y  DIRECTION  OF  CENTER  OF  VIVE  SPECTRUM 
8CDEG=0. 
100     P0  =  90CEC*1. 74532^-? 
XL1=XL*30 
XU1=XU*B0 

call    caulf>T    (xli  »xul » -frol  I  .y.i  cr. 90.  VAR  »T) 
RHsY*0**2*C/S0 
HN=RI1/RI3 

WRITF     (6.50)     BTIDEG.R  I  l.HN.IER 
50     FOP''AT     (•     '.F5.1.T1  7.E)  l.t.T42.F10.<».Tfc2»  13) 
PnnEGrPODEG*?. 

IF  (BOOEG  .LE.  90.)  GO  TO  100 
STOP 
END 
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SUOQorjTT*"7 
THIS    T<? 
1=3 

h=(  ytj-xd/1  :. 

A1=    C.23<i92*'i 

A2  =  ?.<*7862?7 
A3=C.56SP9«? 


•^'TLr^    (XT,1CrJ,t•<:T,Y,T, 


J0,V?  o.T) 


i    ?jpt«    ohpe?    SA'JSS-LFSEWRE    QUAD^AT'JPE 


X  1=-).q0  *17Q  3 


»=. 

'T2  = 
A  =  X 
Y1  = 
Y2  = 

1  =  1 

sn* 

St?* 
B1  = 
B2= 
R3  = 
Rt»  = 
SU 
(F 


5*M 

.5*H 

*2*K2 
+  1 

1  =  3U* 

A-Y1 
A+"2 
J-Y2 


=S0K*fl1*(?"T(31  ,^C,VAP,T)+t"::?(R?rP.0,y5P,T">)  +  A2* 
CT(B3,B0,Vs*,?)+FCT(BU,BQ,Vfi?,T)  )+S3*FCT(  *-  ,R0  ,VAD  ,T) 


A=A+H 
IF(X'J-A)U,U,3 

IF    fSa.l     If3.    i.)    30    ?">    « 
ERP=ABSC  (SC-SOII )  /SU") 

IP    fI.G7.1G)    "0    TO    5 
IF    <FPF-.",1)    ?,c,1 


RETURN 

tt:r=i 
y=su* 

retttr*; 
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t  V      D 


CJf]  D  '"  ~  fj  ""  T  '*"      "T 
'  T  1  y  T  TJ       ">Pp»;        r.PT   'P'^k'T   'j        ^  t»  ri  p  o  v  t  V  I    ""  T  ^  V       <T\  IT 

firrr     ~  t>  »  \t  i^  ^  ^  ^      T  ~,  q  ,f  «  T      ,-. »  0  y  ^ 
T— J* 

i  y  7  ^  ;.*    c*  o  p    wrnsfTV17    V  $  TO  F  S    o  r    / 
rrrv)     1/3,2 

1  1=1 


3  pes     ''»*r)r^ 


in 


V  =  _  "-' 

I1^    ry.LT.    3.°)    m    r 
PHI=1 . 

no  t* n  n  *"■ 

P4I=    .5 

VT  v  (  2  )  =  -  1 .  2  9  3  3  fi  F-  'J 
5"fv  (  3)=-6. 1  S1271E-2 
ft nv  ru)s- 1.287045^-2 

A.  'J  X  (  5  )  =     ^<c"^oc«r.') 
5'l''f7)r     1  .  u  3  3  ?  5  1 T  -  3 


r\  n     1 f     T  =  1    R 

v  -  v  +  v 

?HI=FrIH-Y*ft"H(I  ) 

I?    (I.EQ.    1)    -f!~  =    1.-^tr 
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FUNCTION  FHEAV(3»B0»VAR»T) 
FHEAV=1.595769*(l«*CnSCT*C0 
1  EXP(-<S-B0>**2/(2.*VAl?)) 

RETURN' 
END 


(  B  >  >  >  *  (  1  •  +  C 


1  * 


"  (  3  >  )  )  * 


f unct ion  HEAV(b»BO»VAS*T) 

HEAV  =  1.27324* (l.  +  CGSC  T*C0S(8>)  )*(  1 »  + 

RETURN 

END 


T  *  S I  :\|  t  3  ]   •  l 


rr 


ion 


rorT[  /r  q  r.  y  *  p  < 


1 


FROLL".  39^9U2*  (  1  •+COSCT*COS(B>)  >  *  (  1 
;vp(-(P-pr;)**2/(2»*71s?.)) 

F«D 


.-COS  (T*ST'!  (R)  )) 


p U N  CTT  C  N  8 :  0  T  L  (  9  , *?  0  ,  V  ?. P  ,  T ) 
ROLI=.31<?30^r'*  (L+C^SC^COSCP)  )  )*(  1 
p  £  T  ij  P,  *l 

r,ID 


-COS(T*ST^(B) ) ) 
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Appendix  II.   Directions  For  Use  of  Program 

A.  Discussion. 

This  program  was  intentionally  designed  to  remain  as 
much  as  possible  in  discrete  sections,  each  performing  a 
specific  purpose.   In  this  manner  the  program  can  be  more 
easily  understood  by  others  and  converted  to  use  for  many 
different  sizes  and  types  of  structures. 

An  error  flag  has  been  added  to  the  subprogram  GAULEG. 
If  GAULEG  is  unable  in  a  reasonable  time  to  get  consecutive 
iterations  to  agree  within  one  percent  it  flags  this  fact 
to  the  user  by  indicating  an  error  condition  of  one.   If  the 
error  condition  is  zero  then  GAULAG  has  gotten  the  same  an- 
swer (within  one  percent)  for  two  iterations  with  the 
second  iteration  having  twice  as  many  increments  as  the  first, 

B.  Numerical  parameters. 

To  make  it  easier  to  understand  the  program  flow  the 
leg  spacing  (D) ,  standard  deviation  (SD) ,  and  period  (PER) 
are  simply  assigned  values  in  FORTRAN  statements.   They  could 
just  as  easily  be  read  from  data  cards,  input  from  an  inter- 
active terminal,  or  varied  in  do-loops,  whichever  is  most 
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convenient  for  the  user  on  the  available  system. 

To  adopt  the  program  to  a  square  TLP  of  interest  one 
must  simply  replace  the  assignment  statements  for  leg  spacing 
and  for  the  period  of  interest.   If  many  periods  are  import- 
ant then  a  do- loop  for  period  should  be  used  returning  to 
statement  300  to  ensure  the  value  of  T  is  updated.   If  a 
do-loop  is  used  to  vary  standard  deviation  it  must  return 
to  statement  200  to  revise  the  value  of  the  varience. 

The  center  of  the  wave  spectrum  (3o)  was  varied  in  an 
implied  do-loop  with  statement  100  because  one  is  normally 
interested  in  the  response  to  waves  coming  from  all  directions 
Because  of  the  symmetry  properties  of  a  square  TLP  the  angle 
must  only  be  varied  over  45  degrees  for  heave  and  90  degrees 
for  roll. 

C.   Conversion  to  other  types  of  platforms. 

The  information  concerning  the  platform  geometry  is 
contained  in  the  platform  transfer  functions.   In  the  computer 
program  the  definition  of  the  functions  in  the  function 
subprograms  (ROLL,  FROLL,  HEAV,  and  FHEAV)  entirely  defines 
the  structure.   Note  that  a  factor  of  d2  was  left  in  the  main 
program  rather  than  ROLL  and  FROLL  since  this  factor  will  be 
present  in  moment  transfer  functions  for  any  type  of  structure, 
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Thus  to  change  this  program  to  accomodate  any  shape  TLP 
one  must  simply  perform  derivations  similar  to  those  in  sec- 
tion III-A  and  use  the  results  to  replace  the  function  defin- 
itions in  the  function  subprograms. 

D.   Conversion  to  other  wave  spectra. 

This  program  was  explicitly  written  for  a  normally  dis- 
tributed wave  directional  spectrum  because  it  is  so  easily 
varied.   Simply  by  varying  the  standard  deviation  one  can 
obtain  as  directional  or  uniform  a  spreading  as  desired. 
However,  if  another  spectrum  is  preferred  this  program  can 
be  adapted. 

The  most  obvious  place  to  change  is  in  the  definition 
of  FROLL  and  FHEAV.   There  the  new  spectrum  would  be  substi- 
tuted for  the  exponential  term.   In  the  main  program  the 
correction  factor  C/SD  is  applied  after  calling  GAULEG  using 
FROLL  or  FHEAV.   This  would  be  changed  or  eliminated.   Since 
C,  the  factor  for  loss  of  tails,  is  no  longer  necessary  the 
entire  subroutine  ERF  could  then  be  eliminated.   New  correc- 
tion factors  for  the  new  spectrum  would  then  be  introduced. 

It  is  also  possible  to  normalize  the  transfer  function 
using  something  other  than  uniform  spreading.   For  this  ROLL 
and  HEAV  only  must  be  revised.   If  correction  factors  are 
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necessary  they  are  best  introduced  in  the  main  program  after 
calling  GAULEG  for  the  integration. 
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Appendix  III.   Gauss -Legendre  Quadrature 

A.   Legendre  Polynomials 

1.   Definition.   Define  a  set  of  polynomials  Pn(x)  of 
degree  n,  where  n  is  a  non- negative  integer,  such  that 


.{  pm<*>  Vx)  dx  =  cn  6mn 


and       P„(D  s  1. 
n 


The  first  is  an  orthogonality  condition  and  the  second 
is  necessary  to  fix  cn.   Together  they  define  a  unique  set  of 
polynomials  known  as  the  Legendre  polynomials,  the  first  few 
of  which  are 

P0(x)  =  1 

P1(X)   ■  X 

P2(x)  -  J5<3X2-1) 

P3(x)  «  J5(5xs-3x)  . 

In  general  they  can  be  found  by 

Pn(x)  =  -i—  3IL  (x2-l)n. 
n      2   n!  dxn 

Note  that  P  (x)  has  n  zeroes  in  the  interval  from  minus  one 
n 

to  one  and  is  even  for  n  even  and  odd  for  n  odd. 


67 


2.   Polynomial  approximation.   Any  function  piecewise 
continuous  in  the  interval  from  minus  one  to  one  can  be 
expanded  in  terms  of  an  infinite  series  of  Legendre  polynomials 
It  can  be  approximated  by  a  finite  series  of  Legendre  poly- 
nomials. 

y(x)  *mIoCm  Pm<X)  SnJ0Cra  Pm(X> 


B.   Integral  approximation 

A  simple  means  to  evaluate  the  integral  in  the  region 
minus  one  to  one  is  as  follows: 

1  1   OS 

/  y(x)  dx  »  /  I  c-  pm<x)  dx 

-l  -10 

Multiply  both  sides  by  1  *  PQ(x): 

}  y(x)  dx  »  I  c   }  Pm(x)  P  (x)  dx 
c0  /  PQ(x)  PQ(x)  dx 

/  y(x)  dx  =»  2  c_ 

-i  ° 

Thus  the  problem  has  now  become  one  of  finding  c  ,    the 
coefficient  of  P  (x) .   One  possible  method  will  be  shown 
using  the  example  where  n=2. 

We  desire  to  find  coefficients  and  arguments  such  that 
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_|    y(x)    dx   =    1   C0    =   Afl    y(x0)    +   A^    yU^ 


For   any    two  points.    x0    and   xx ,    we   expect 

Noting   that  P.  (x)    ■    1   and  P.  (x)    =   x   and   choosing   x      and 

0  1  a    0 

x   to  be  the  zeroes  of  P  (x)  this  becomes 


Co  +  Ci  Xo  =  ir(xo) 


C   +  C   x   =  y(x  ) 


Where  x   and  x   are  ±  (1//3)  .   By  adding  these  two 
equations  and  noting  that  xfl  ■  -xJ  we  have 


2  CQ  =  y(-l//3)  +  y(l//3) 


which  is  our  desired  answer. 

To  get  a  larger  degree  of  accuracy  a  larger  n  is  required, 
This  method  for  determining  2  C0  remains  effective.   For  the 
case  of  n  =  5  the  solution  is  (to  seven  significant  digits) 

i 

Jy(x)  dx  ■  (0.568  888  9)  {y(0)> 

+  (0.478  628  7)  {y(.538  469  3)+y(-.538  469  3)} 

+  (0.236  926  9)  {y(.906  179  8)+y(-.906  179  8)}  . 
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However,  one  is  not  always  interested  in  the  interval 
from  minus  one  to  one.   A  transform  must  then  be  used  to 
transfer  the  region  of  interest  to  that  interval.   For  addi- 
tional accuracy  without  the  need  for  a  larger  order  quadra- 
ture the  total  interval  can  be  divided  into  subintervals. 
Each  subinterval  can  be  transformed  seperately  and  the 
result  summed.   Thus  each  subinterval  will  be  approximated 
by  a  fifth  order  polynomial.   This  is  the  method  employed 
in  the  subroutine  GAULEG. 
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